Optical Generation of Vortices in trapped Bose-Einstein condensates 
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We demonstrate numerically the efficient generation of vortices in Bose-Einstein condensates 
(BEC) by using a "phase imprinting" method. The method consist of passing a far off resonant 
laser pulse through an absorption plate with azimuthally dependent absorption coefficient, imaging 
the laser beam onto a BEC, and thus creating the corresponding non-dissipative Stark shift potential 
and condensate phase shift. In our calculations we take into account experimental imperfections. We 
also propose an interference method to detect vortices by coherently pushing part of the condensate 
• using optically induced Bragg scattering. 
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One of the remaining challenges of the physics of 
, trapped Bose-Einstein condensates (BEC) con- 
^ ' cerns the demonstration of their superfLuid behavior. Su- 
perfluidity is inevitably related to the existence of vor- 
J2 ■ tices and persistent currents in BEC, which so far have 
not been observed experimentally, despite serious efforts 
jH . Two aspects of the vortex problem have been studied 
, intensively: i) in rotating traps vortices appear in a natu- 
Q ■ ral way as thermodynamic ground states with quantized 
^ O ^ [ angular momentum . Stability and other properties of 
vortices in rotating traps have been throughly discussed 
I '. in Ref. (7|-|^; ii) in stationary traps creation of vortices 
K*' ' (or related dark solitons in ID) requires the use of dy- 
namical means, and an independent stability analysis. 
^1 ■ Several methods were proposed to generate vortices 
J — , in non-rotating traps: stirring the condensate using a 
■ blue detuned laser [fo 11 1, or several laser beams 



0^ ' adiabatic passage |l^ or Raman transitions in bi- 
condensate systems. Such vortices are typically not sta- 
ble, and can exhibit dynamical or energetic instability. In 
the first case vortices decay rapidly, in the second the vor- 
^ ] tices are stable within the framework of the mean field 
■ theory, and their corresponding decay requires to take 
^ I into account interactions of the BEC with the thermal 
O ■ cloud. In the latter case the vortex dynamics is expected 
to be sufficiently slow, and thus experimentally accessible 

k>( _ In this Communication we propose and investigate yet 
; I ■ another procedure of vortex generation using "phase im- 
[ printing". This method consists of i) passing a far off 
resonant laser pulse through an absorption plate whose 
absorption coefficient depends on the rotation angle ip 
around the propagation axis and ii) creating the corre- 
sponding Stark shift potential inside a BEC by imaging 
the laser pulse onto the condensate which leads to a </? 
dependent phase shift in the condensate wave function. 
This method is very efficient and robust, and allows for 
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engineering of a variety of excited states of BEC con- 
taining vortices. In the ideal case the method allows to 
generate genuine vortices with integer angular momenta. 
In the presence of imperfections, typically more complex 
vortex patterns are generated. 

We suggest an interference based method of vortex de- 
tection, in analogy to methods used in nonlinear optics 
| p6[ . In the context of matter waves similar methods 
were proposed in references |l^Jl^. Our idea is to com- 
bine this method with the recently developed techniques 
of Bragg diffraction for BEC manipulation [|l9| . 

Before we turn to details, we should stress that the 
dynamical generation of vortices differs from the case of 
rotating traps, in which a pure vortex state with angu- 
lar momentum = 1 (in units of h) is selected in the 
process of reaching the equilibrium. In our case, genera- 
tion of pure vortices a requires fine tuning of parameters 
which is hard to achieve in experiments. Our method is 
suitable for creation of generic states with vorticity , 
i.e. states with several vortex lines, around which the 
circulation of velocity does not vanish |Q . 

Dynamics of BEC, and thus the process of creation 
and evolution of vorticity at zero temperature is well de- 
scribed by the time-dependent Gross-Pitaevskii equation 
for the wave function ip{r, t): 



where Vt(r) = M{uj1x'^ + ujy-y'^ + cj^z^)/2 is the external 
trap potential which we assume to be harmonic, M is the 
mass of the atom, and LOx,^y, are the trap frequencies. 
The nonlinear term VAri(r,t) — g|'i/;(r,t)p describes the 
mean field two-body repulsive interaction whose strength 
g is related to the scattering length a by g = ^■KNfP'a/M , 
where N is the total number of condensed particles. The 
term V;(r,t) describes an effective potential created by 
an external laser beam impinging on the condensate af- 
ter passing through a plate with appropriately modulated 
absorption coefficient. In the following we will study the 
two dimensional version of the Eq. (|l|), in which we re- 
place g hy g/D, where D is the characteristic depth cho- 
sen such that the ground state chemical potentials of the 
3D and 2D systems are equal. 

In this Communication we show that a short pulse of 
light with a typical duration of the order of fractions of 
microseconds with properly modulated intensity profile 
creates vorticity in a Bose-Einstein condensate initially 
in its ground state. If the incident light is detuned far 
from the atomic transition frequency its main effect on 
the atoms is to induce a Stark shift of the internal energy 
levels. As the intensity of light depends on the position, 
the Stark shift will also be position dependent as will be 
the phase of the condensate. 

The main feature characterizing a vortex is related to 
the particular behavior of the phase of the wave-function 
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at the vortex line: the phase "winds up" around this hne 
i.e. it changes by an integer multiple m of 27r on a path 
surrounding the vortex. Index m is the vortex charge. 
The light beam, before impinging on the atomic system, 
is shaped by an absorption plate whose absorption coef- 
ficients varies linearly with rotation angle (f around the 
plate axis (see Fig. 1). As a result, Vi{r, t) depends on the 
distance from the propagation axis p, and the azimuthal 
angle (p. 

In the ideal case this absorption plate causes a real 
jump of the potential at, say ip = 0. In this case we can 
model the potential for < t < T by 



being zerofor other times. Here T is the (square) pulse 
duration, whereas / denotes the characteristic Stark shift 
(proportional to the laser intensity); L is the charac- 
teristic length scale, on which the absorption profile 
is smoothed in our calculations in the vicinity of the 
propagation-rotation axis. The use of more realistic tem- 
poral pulse shapes does not change the results signifi- 
cantly. In the following we will use square pulses, since 
it allows to control the pulse area in a simpler manner. 

The characteristic time and length scales for the BEC 
that we consider are ms, and /zm, respectively Note, 
that for this parameter range T is short (~ fis), and 
III much larger then other energy scales, the dynamical 
effect of Vi (jO, tp) corresponds approximately to a "phase 
imprinting" . Since during the interaction with the laser, 
all other terms in Eq. (||) can be neglected, the wave 
function after switching off the pulse becomes 



For the ideal case, choosing IT = to, we get the desired 
phase dependence characteristic for a pure vortex state. 
Unfortunately, the above picture is oversimplified, be- 
cause the wave function (||) does not necessarily vanish 
at the vortex line {p = 0), which signifies infinite kinetic 
energy and its unphysical character. One has to solve 
the full dynamics of the system for < i < T, taking es- 
pecially into account the kinetic energy term. Moreover, 
realistic absorption plates cannot generate the singular- 
ity in the i^j-dependence of Vi . For this reason we smooth 
the (/^-dependence of the intensity / on the scale of S radi- 
ans (see Fig. 1). "Phase imprinting" with this smoothed 
intensity distribution cannot generate pure vortex states 
— it can, however, as we shall see below, generate in a 
controlled way states with predetermined vorticity. 

The Gross-Pitaevskii equation was solved using the 
split operator method in 2D {x and y plane). The sim- 
ulations were divided into two stages: the initial exci- 
tations stage of duration T ~ 0.16fis, with about 1000 
time steps, and the second stage with a characteristic 
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V'(r, T) = exp{-iTVi{p, <p)/h)i:{r, 0). 



(3) 
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time scale of ms, and about 1000 steps/ms. We have as- 
sumed that initially a condensate containing TV = 100000 
Rubidium atoms (with a = 5.8nm) was trapped in a 
disc-shaped trap of frequencies — cijj, = 27r x 30Hz, 
bJz — 2n X 300Hz. In this case the condensate radius 
was about 15/im. We have used / = 6 x lO^Hz, and ad- 
justed T to obtain the desired values of the pulse area 
{IT ~ 1,2, . . .). The /9-dependence of the potential Vi 
was smoothed with L ~ 1.7/im, and the (^-dependence 
with 5 ~ 0.15 radians. The smallest possible value of 
5 ~ 0.04 radians was determined by the spatial grid size 
which typically was 512x512 points in a 40/zm x 40/im 
box. Additionally, we studied the case when the laser 
was focused slightly off the trap symmetry axis. During 
the dynamics, after switching off the laser, the following 
quantities are conserved: wave function norm, mean en- 
ergy, and mean z-component of the angular momentum. 
They were monitored in order to control the accuracy of 
numerics. The first two quantities were constant within 
the accuracy of the method, the third one in some cases 
exhibits slow variations of the order of a percent. 

One of the most important questions concerning the in- 
vestigation of vortices is an efficient method for their de- 
tection. Experimentally monitoring density profiles with 
the necessary resolution is difficult, since the vortex core 
is very small. The best way is to monitor the phase of 
the wave function in an interference measurement. Such 
interference measurements are routinely done in nonlin- 
ear optics [|l^. In the context of BEG, they were first 
proposed by Bolda and Walls jl7| , who considered the in- 
terference of two condensates moving towards each other. 
If both condensates are in the ground state (no vortices), 
one expects interference fringes as those observed by Ket- 
terle et al. |2^. In the case of interference of one con- 
densate in the ground state with the second one in the 
m — 1 vortex state, a fork-like dislocation in the inter- 
ference pattern appears. The distance between the inter- 
ference fringes is determined by the relative velocity of 
the condensates, which can be controlled experimentally. 
This is a very efficient and clear method of vorticity de- 
tection. It requires, however, the use of two independent 
condensates. 

We propose here to combine the interference method 
with the recently developed Bragg diffraction technique 
, which has been successfully used in a four wave mix- 
ing experiment ||2^. The idea is to transfer part of the 
atoms coherently to another momentum state using one, 
or several stimulated two-photon Raman scattering pro- 
cesses. The procedure is the following: i) first we create 
the vortex, or vorticity state in the trap; ii) we open the 
trap and let the condensate expand; iii) when the den- 
sity is reduced to values for which nonlinear interactions 
are negligible, we apply the Bragg pulses. Part of the 
wave function attains a phase factor that signifies the 
fact that the corresponding momentum was transferred 
to part of the atomic sample pj]. The resulting wave 
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function is the superposition of two vortex (or vorticity) 
states moving apart from each other, with a velocity that 
can be easily controlled by choice of the angle between 
the Bragg beams. In the following we will use velocities 
of the order of 1 mm/s, which allow for efficient detec- 
tion after 4-5 ms when the vortices are about 5/xm apart; 
iv) detection consists of optical imaging that is accom- 
plished within a few /xs. The interference patterns have 
a characteristic length scale of a few fim. 

We shall first discuss our results for the ideal case: 
5 < 0.04 radians, IT = 1, and the laser focused exactly 
along the trap symmetry axis. In this case during the first 
stage of the dynamics angular momentum is transferred 
to the system. The mean of Lz grows monotonically from 
to about 1 (in units of h), and stays close to 1 hereafter. 
The dispersion SLz grows within the pulse duration to 
about 5, returns to about for t = T. In the second stage 
of the dynamics, a hole in the density profile appears, and 
the density acquires a stationary state within a fraction 
of a ms. This time scale corresponds to the inverse of the 
characteristic frequency, determined by the mean field 
energy in the condensate, g/hV, where V is the volume 
of the condensate {V ~ 7rl5 x 15 x 4/im^). It is the time 
required by the second sound to travel over a distance of 
the order of the healing length. The density distribution 
hereafter remains constant on the scale of 10 — 20ms. 
The dispersion of the angular momentum, SL^, exhibits 
a stepwise growth to about SL^ ~ 5 with fast growth 
periods of duration of about 0.5 ms, followed by slow 
growth stages of duration of a few ms. We observed a 
slight movement of the vortex off the center of the trap. 

In Fig. 2a we present a typical simulation result ob- 
tained a few ms after the creation of a vortex using the 
detection method of based on the interference of a 
condensate in the vortex state with a condensate in the 
plane wave state. The characteristic fork-like pattern 
reflects the fact that the phase winds up by 27r as one 
circulates around the vortex line. Fig. 2b shows the 
simulation of the interference pattern for the method of 
detection proposed by us. A few milliseconds after ap- 
plying the Bragg pulses we can clearly see a double fork 
structure, as one should expect for this case. The forks 
are oriented in opposite directions, because the conden- 
sates have the same helicity, but opposite velocities. 

It is interesting to compare this ideal result with more 
realistic simulations obtained for S = 0.15, keeping the 
laser focused at the center of the trap and the pulse area 
IT =1. In this case one excites a vorticity state which 
is not pure. The mean after the laser pulse is close to 
zero, as it should be in the absence of the phase discon- 
tinuity. The dispersion SL^ behaves similarly as before, 
although it grows a little faster. The density profile ex- 
hibits now the hole in the center but also a number of 
deep minima, or holes at iy9 = in the outer regions of 
the condensate. While the density distribution is difhcult 
to interpret, the interference method allows to interpret 
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them very clearly. Application of the plane wave method 
indicates two vortices of opposite winding number (topo- 
logical charge), and perhaps two more close to the edge 
of the condensate (Fig. 3a). This basic structure is dou- 
bled, but nevertheless clearly recognizable using interfer- 
ence combined with the Bragg scattering (Fig. 3b). 

Finally, in Fig. 4 we present numerical results that 
mimic other possible experimental imperfection: offset of 
the laser focus by 1.7/im off the center of the trap (ver- 
tical in the figure plane), combined with the smoothed 
</?-dependence of absorption profile {S = 0.15). In this 
case the behavior of and SL^ is similar to the case 
presented in Fig. 3. The interference patterns, obvi- 
ously, are more noisy and irregular. Still one can clearly 
recognize fork-like dislocations corresponding to vortices 
using the plane wave interference method (Fig. 4a), and 
the method proposed here (Fig. 4b). 

In our numerical simulations we have also studied sit- 
uations when the area of the phase imprinting pulse is 
large, IT > 1. Typically a number of singly charged 
vortices aligned in geometric patterns similar to those 
obtained in rotating traps Q are created. 

Summarizing, we have proposed a method of vortex 
generation in trapped Bose-Einstein condensates. The 
method employs optical potentials induced by passing a 
laser beam through an absorption plate with an absorp- 
tion coefficient that depends on the azimuthal angle ip. 
We have also proposed a method of detecting vortices and 
vorticity by combining the known interference method 
with the recently developed Bragg scattering techniques 
to coherently transfer atoms into a selected momentum 
state. 
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Krolikowski, and K. Rzazewski. This work is supported 
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FIG. 1. Absorption plate. L is the radial extension on 
which the absorption profile was smoothed in our simulations, 
whereas 5 is the angular extension of the smoothing. 




FIG. 2. Interference pattern of the vortex state (5 ms after 
phase imprinting): (a) superimposed with a plane wave of 
k — 1.8 /im~^ moving vertically in the plane of the figure; 
(b) 5 ms after applying Bragg pulses transferring momentum 
to a part of the condesate (along vertical axis in the plane 
of the figure). The ideal case of equal splitting of the whole 
condensate has been assumed. The incident pulse was focused 
at the trap center and had a sharp step in the intensity profile 
{S = 0.04). 




FIG. 3. Interference pattern of the vorticity state (6 ms 
after the phase imprinting): (a) after superimposing with a 
plane wave as in Fig. la; (b) after applying the Bragg pulses 
as in Fig. lb. The incident pulse was focused at the trap 
center and its intensity profile was smoothed over an angle 
S = 0.15. 




FIG. 4. Interference pattern of the vorticity state (8.5 ms 
after phase imprinting) : (a) after superimposing with a plane 
wave as in Fig. la; (b) after applying the Bragg pulses as in 
Fig. lb. The incident pulse was focused 1.7/xm off the trap 
center and its intensity profile was smoothed over an angle 
S = 0.15. 
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